% Written by: Orhan Torul
% Istanbul, Turkey, May 2017
function Gvec = sgump(varmp,fy1,fy,fx1,fx)

global ga ha;

a11_g=varmp(1);
a12_g=varmp(2);
a21_g=varmp(3);
a22_g=varmp(4);
a31_g=varmp(5);
a32_g=varmp(6);
a41_g=varmp(7);
a42_g=varmp(8);
a51_g=varmp(9);
a52_g=varmp(10);
a61_g=varmp(11);
a62_g=varmp(12);
a71_g=varmp(13);
a72_g=varmp(14);
a81_g=varmp(15);
a82_g=varmp(16);
a91_g=varmp(17);
a92_g=varmp(18);
a11_h=varmp(19);
a12_h=varmp(20);
a21_h=varmp(21);
a22_h=varmp(22);
f_y_1=fy1;
f_y=fy;
f_x_1=fx1;
f_x=fx;
ga=[a11_g a12_g; a21_g a22_g; a31_g a32_g; a41_g a42_g; a51_g a52_g; a61_g a62_g; a71_g a72_g; a81_g a82_g; a91_g a92_g];
ha=[a11_h a12_h; a21_h a22_h];

Gara =f_y_1*ga*ha + f_y*ga + f_x_1*ha + f_x;
Gvec=Gara(:);
g_matrix_mp=ga;
h_matrix_mp=ha;
